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ABSTRACT 



Several recent studies have reported the detection of an anomalous color spread along the red giant branch (RGB) of some globular 
clusters (GC) that appears only when color indices including a near ultraviolet band (such as Johnson U or Stromgren u) are consid- 
ered. This anomalous spread in color indexes such as U-B or c y has been shown to correlate with variations in the abundances of light 
elements such as C, N, O, Na, etc., which, in turn, are generally believed to be associated with subsequent star formation episodes that 
occurred in the earliest few 10 s yr of the cluster's life. Here we use publicly available u, g, r Sloan Digital Sky Survey photometry to 
search for anomalous u-g spreads in the RGBs of nine Galactic GCs. In seven of them (M 2, M 3, M 5, M 13, M 15, M 92 and M 53), 
we find evidence of a statistically significant spread in the u-g color, not seen in g-r and not accounted for by observational effects. In 
the case of M 5, we demonstrate that the observed u-g color spread correlates with the observed abundances of Na, the redder stars 
being richer in Na than the bluer ones. In all the seven clusters displaying a significant u-g color spread, we find that the stars on the red 
and blue sides of the RGB, in (g, u-g) color magnitude diagrams, have significantly different radial distributions. In particular, the red 
stars (generally identified with the second generation of cluster stars, in the current scenario) are always more centrally concentrated 
than blue stars (generally identified with the first generation) over the range sampled by the data (0.5r/, $ r S 5r h ), in qualitative 
agreement with the predictions of some recent models of the formation and chemical evolution of GCs. Our results suggest that the 
difference in the radial distribution between first and second generation stars may be a general characteristic of GCs. 
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1. Introduction 



stars: red giant branch -GCs: individual (M 2, M 15, M 13 , M 5, M 3, M 53, M 92, NGC 2419, 



Globular clusters (GCs) have long been considered as examples 
of the theoretical concept of a simple stellar population, i.e., 
a population o f stars that is strict l y coe val and chemically 
homogeneous dRenzini & Buzzonil Tl986h . This traditional 
paradigm remains valid in some cases, but recent photometric 
and spectroscopic observational results strongly indicate that 
most GCs have been sites of two or even more star-formation 
episodes, p roducing a peculiar chemical (s elf-)enrichment 
pattern (see iPiottol l2009t ICarretta et al.ll2010dl and references 
therein, for review and discussion). 

In general, GCs are largely homogeneous if we consider 
Fe-peak elements but exhibit a significant spread in the 
abundance of lighter elements, with strong anti-correlations 
between, for example, the abundances of Na and O, or Mg 
and Al, as w ell as bimodal distributions of CH and CN line 
strengths (see ICarretta et aT1l2010dt iMartell & Smithll2009l and 
references th erein). While these anomalies have been known for 
decades (see lKraf3ll994 for a review of early results), it has 
only recently been establish ed that (a) they can also be traced 
in un-evolved stars ( Gratton et alj 1200 U iRamfrez & Cohenl 
120021: ICarretta et alJ 120041) . hence cannot be due to mixing 
processes, and (b) they seem to be a characteristic feature 
of GCs ( ICarretta et alj 1201 Odl) . Interestingly enough, these 
chemical inhomogeneities are not confined only to Galactic 



GCs dColucci et alll2009l:lMucciareili et al.lr2009l:lJohnson et all 
2006: lLetarteetalJ f2006. and references therein). 



The scenario generally invoked to explain the above phe- 
nomena, foresees two subsequent generations of stars in 
the first few hundred Myr of cluster life, with ejecta from 
the massive stars of the first generation enriching the in- 
tra cluster medium (ICM) before the formation of the sec- 
ond on e. Intermediate-mass asymp totic giant branch stars (1M- 
AGBs; Ventura & D' Antona 2008) or fast rotating massive stars 
(FRMSs ; lDecressin et alj|2007l) have been proposed as the ICM 
polluters that most likely explain the observed patterns (see 
Renzini 2008, for a discussion of the merits of the two models). 

In this context, several results imply that the spread in light- 
element abundances of GCs can be traced by photometric in- 
dices with near ultraviolet (wide or intermediate) passbands, en- 
compa ssing the wavelen gth range 3000A55 A S= 4000A. For in- 
stance, Yong et al. (2008) demonstrated that the c v index (includ- 
ing the near-ul traviolet Stromgren u passband, akin to the c\ in- 
dex used bv iGrundahl & Brilevl [20011: IGrundahl et alj|2002l to 
trace NFQ) very clearly traces the differences in N abundances 
among the stars of NGC 6752, over a range of ~ 2.0 dex in 
[N/Fe]. As the abundance of N is found to correlate (or anti- 
correlate) very well with other light elements exhibiting abun- 



1 In particular, c y = C\—(b—y) and C\ 
(2008), and references therein. 



(u-v)-(v-b), see Yo ng et all 
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dance spreads in GCs (O, Na, Mg, Al), differences in these 
elements can also be co r related with c y variations, as demon- 
strated by ICarretta et al.l d2009bl) . These authors used their Na 
and O abundance determinations for NGC 6752 to show that 
stars with different abundances of these elements lie on differ- 
ent sides of the cluster RGB in the (V,c y ) plane, and concluded 
that Na-poor (Na-r ich) stars are also N-poor (N-rich; see also 
iMilone et all 1201 Pi for further results on NGC 6752 and dis- 
cussion). Moreover, the well-known anti-correlation between Na 
and O implies tha t Na-rich stars are also O-poor and vice versa. 
Yong et all d2008[) noted that available Strom gren photometry of 
Galactic GCs dGrundahl & Andersen| [T999) suggests that a sig- 
nificant spread in N abundance may be a general property of 
GCs. In particular, their Fig. 14 shows remarkable spreads in c y , 
at any magnitude along the RGB of NGC 288, NGC 362, M 5, 
M 3, M 13, NGC 6752, M 15, M 92, and NGC 6397. 

Mari no et al.l (2008) were able to demonstrate for M 4 that 
a bimodal Na-O distribution, correlated with a bimodal distribu- 
tion in CN strength, was also clearly associated with a bimodal 
spread in the color of RGB stars in the U vs. (U-B) CMD, not 
seen with other color indices. In this case, as in all cases in which 
the U-B color spread in the RGB has been correlated with O, Na 
abundances, the Na-rich/CN-strong stars lie on the red side of 
the RGB, while Na-poor/CN- weak stars populat e a bluer por- 
tion of the branch. In line wifh lYong et all (12008). M arino et all 

(2008) attribute the spread in (U-B) color to strength variations 
of several NH and CN bands included within the U passband, in 
particular the NH band around 3360A, and the CN bands around 
3590A, 3883A and 4215A. They performed synthetic U, B pho- 
tometry of theoretical spectra with different strengths of the NH 
and CN features, mimicking the observed spectra of M 4 stars, 
finding that the resulting variation in the U-B color goes in the 
right direction but has an amplitude four times smaller than the 
observed color spread. 

iKravtsov et al.1 (l2010allbh used U band photometry to 
identify multiple populations within the clusters NGC 3201 
and NGC 1261. The spread in U-B color among RGB stars 
of NGC 3201 was sh own to c orrelate wi t h Na abundance by 
ICarretta et ail d2010dl) . Finally, lHan et al.l d2009l) showed that 
the RGB of NGC 1851, which is narrow and well defined 
in optical CMDs not including U photometry, is cl early split 
into tw o parallel branches in the U vs. U-I colour. lHan et all 

(2009) suggest that the splitting is caused by a combination 
of effects due to variations in the abundances of not only of 
C,N,0 but also of heavier elements (Ca, Si, Ti, and Fe) and 
helium. The presence of a small (but real) spread in iron and 
Q'-elements abundan ce in this clu s ter wa s recently confirmed 
spectroscopically by ICarretta et al.l (1201 Pel see also their Fig. 4, 
where they show that the cluster stars segregate along the RGB 
according to their Na abundance as in the other cases described 
above). Additional results and a more general discussion about 
the possible role of heavy elements variati ons in producin g the 
RGB color spreads in GCs can be found in lLee et all d2009l) and 
ICarretta etaf]d2010bl) . 

While the details of the process remain still unclear, it seems 
reasonably well established that (at least) light-element abun- 
dance spreads in GCs are correlated with spreads in the colors 
of RGB stars if near ultraviolet filters are adopted. This opens 
an interesting window in investigating the origin of the chem- 
ical evolution of GCs as accurate U photometry for very large 
sample of GC stars can be obtained in a much easier way and 
for much more distant clusters than the mid-to-high-resolution 



spectroscopy needed to obtain direct chemical abundance esti- 
mates (which, however, provides more information). In partic- 
ular, wide-field photometry capable of discriminating between 
N-rich and N-poor stars would provide the large samples and 
the wide radial coverage that are needed to compare the radial 
distributions of the two groups; this may provide a very interest- 
ing insight into the process of GC formation and early chemi- 
cal enrichmen t (see lD'Ercole et al]|2008t [becressin et al.| [2007; 
lRenzinill2008l) . 

Prompted by these considerations we decided to test whether 
additional observational evidence of color spreads in the RGB of 
GCs could be derived from the publicly ava ilable Sloan D igital 
Sky Survey photometry (SDSS; see lAbazaiian et aLll2009l and 
references therein), as it provides multicolor photometry includ- 
ing a u passbancQ and incorporates several Galactic GCs. In 
the following, w e use the accurate ugriz SDSS photometry by 
lAn et alJ d2008l) of nine GCs selected from their sample to show 
that: (a) the u-g color index indeed correlates with spectroscopic 
Na,0 abundances and that the associated spread in color can 
be detected in SDSS data, at least in clusters of intermediate to 
high metallicity; (b) a statistically significant intrinsic u-g spread 
among RGB stars is detected in seven of the nine clusters consid- 
ered in the present study; and, finally, (c) in these seven clusters 
the radial distribution of stars lying on the blue side of the RGB 
(in u-g) differs significantly from that of stars lying on the red 
side, the latter being more centrally concentrated than the for- 
mer. 

In Sect. |2] we briefly describe the adopted photometric data 
set and present evidence of a correlation between light-element 
abundances and u-g spread along the RGB for the cluster 
M 5. In Sect. [3] we analyze the g,g-r and g,u-g CMDs of the 
considered clusters in search of u-g color spread in the RGB 
that cannot be accounted for by observational effects. We finally 
discuss our results in Sect. [4] Some preliminary results from 
this study were presented by M. Bellazzini at the meeting "The 
Giant Branches" held in Leiden in May 20090- 



2. The data set 

lAn et al.l d2008l An08 hereafter) reanalyzed SDSS images of 
the GCs (and open clusters) included in the survey using 
the D AOPHOT/ALLFRAME suite of programs dStetsonlll987l 
1994). These programs are known to perform more effectively 
than the survey photometric pipeline in the high-crowding con- 
ditions typical of dense star clusters. 

We limited our analysis to the most favorable cases, exclud- 
ing clusters with \b\ < 20° - to avoid strong contamination from 
Galactic field stars - and/or with E(B - V) > 0.15 - i.e., clusters 
with relatively high extinction. Moreover, we decided to include 
in our sample only clusters with more than 100 candidate RGB 
stars between the horizontal branch (HB) and two magnitudes 
below this level, to ensure a solid statistic basis to our analysis. 
According to these criteria, among the seventeen clusters con- 
sidered by An08 we selected the nine GCs listed in Tab.Q] 

Only stars with valid magnitude estimates in both g and 
r have been retained for subsequent analysis. When An08 

2 The SDSS u filter has A eff =352\A and FWHM=555A. 
For details about the photometric filters quoted in this 
paper see the Asiago Database of Photo metric Systems 
http : //ulisse .pd. astro. it/Astro/ADPS/ (Moro & Munari 
12050) . 

3 www . lorentzcenter . nl/lc/web/2009/324/Friday/Bellazini . ppt 
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Fig. 1. g,g-r and g,u-g CMDs for M 5, M 15. The stars for which 
spectro scopic abundances are available from from lCarretta et all 
(2009ab) are marked according to their Na abundances: Na-poor 
stars are marked as blue open circles, while Na-rich stars are 
marked with red open triangles. Stars plotted as heavy dark dots 
are those selected as candidate RGB stars from the g,g-rCMD. 



provided more that one catalog per cluster (i.e., in cases of 
clusters imaged in d ifferent overlapping SDSS "plates", see 
Abazaiian et al. 2009, An08, and references therein), we merged 
the catalogs into a single one including all the cluster stars listed 
by An08 without duplications. The final catalogs cover fields in- 
cluding the vast majority of cluster stars but does not reach the 
tidal radius (except for NGC 2419, see below). 

In some cases, the available photometry was highly incom- 
plete and/or was characterized by large uncertainties in the 
innermost regions of the cluster, because of the high degree 
of crowding. In most cases the final catalog also covered a 
field of view including regions where the contamination from 
Galactic stars was not negligible. To avoid the inclusion of clus- 
ter stars with poor-quality photometry and contamination from 



fore/background Galactic stars, we limited the analysis to a ra- 
dial corona between the r,„ and r oul values specified in Tab. Q] 
In this table, r- m and r out are also expressed in units of core, 
half-light, and tidal radii (r c , r/,, and r t , respectively), to pro- 
vide a more objective idea of the actual radial sampling and of 
the cluster-to-cluster differences in the radial sampling. We note 
that: (a) in all the considered cases we sample the clusters out 
to > 4 rh, and in most cases out to 6 - 10 r/,; (b) in the cases 
in which we excluded the innermost 1' from the analysis, the 
innermost regions between ~ 1 - 3 r c were lost, the complete- 
ness in these inner parts possibly being remarkably low also in 
clusters for which we retained the central regiotQ and (c) the 
considered clusters span a wid e range of ce ntral density (more 
than 4 orders of magnitude, see iHarrisI 19961) . hence the different 
radial ranges also sample regions in widely different dynamical 
conditions, depending on the considered cluster. 

The reddening for the GCs are taken from lHarrii dl996l) . The 
coefficients of the adopte d extinction laws (A ^/Ay) are taken 
from the computations by Girardi et al] d2004l) for cool giants 
(r c// =4000, log g=2.00 and [M/H]=-2), and assuming A v = 
3AE(B - V). In particular: A u = 4ME(B - V), A g = 3.64£(fl - 
V), and A r = 2.1\E(B - V), and E(u - g) = 1.20£(B - V) and 
E(g - r) = Q.93E(B - V). 



2.1. Na abundances and u-g colors in RGBs 

For two of the selected clusters, M 5 and M 15, we found a 
significant sample of stars for which there are spectroscopic Na 
abundances over a remarkably wid e luminosity range along the 
RGB, from lCarretta et all (l2009albh . Following the criteria used 
by these authors, we divide the RGB stars between a candi- 
date first generation and a candidate second generation (P and 
I+E components, respectively, adopting their nomenclature^) 
by adopting a threshold in sodium abundance that depends on 
the cluster metallicity. In particular, stars having [Na/Fe] m i„ < 
[Na/Fe] < [Na/Fe] m j„ + 0.3 are assigned to the Na-poor P com- 
ponent, while stars having [Na/Fe] > [Na /Fe] m i„ + 0.3 are as- 
signed to the Na-rich I+E component (see Carre tta et alJ l2009b. 
for further details). The resulting thresholds are [Na/Fe] = +0. 10 
for M 5 and [Na/Fe] = +0.22 for M 15; we stress, however, that 
the results presented below are not particularly sensitive to the 
adopted threshold. 

The upper panel of Fig. Q] shows that, in the case of M 5, 
the Na-poor and Na-rich stars, that are tightly aligned along 
the narrow cluster RGB in the g,g-r CMD, are clearly separated 
into two parallel sequences in the much broader giant branch 
seen in the g,u-g diagram, the Na-rich stars appearing system- 
atically redder than Na-poor ones (at least below the HB l evel), 
a beha vior strictly analogous to that o bserved by iMarino et alj 
(120081) in M 4, bv lMilone et all d2010h in NGC 6752, using U 
photometry, and bv lCarretta et al.1 (l2009bl) . again in NGC 6752, 
but using Stromgren u photometry. 

This result clearly indicates that SDSS u photometry can 
also be used to trace the UV color spread correlated with the 
light-element abundance spread described in Sect. [TJ above. We 



4 This should not affect the results presented below, as in all the cases 
we compare RGB stars having the same distribution in magnitude, see 



Sect.[3TTI 

5 P stands for the first - Primordial - generation (Na-poor stars), while 
Na-rich stars are divided into I (Intermediate) and E (Extreme) subse- 
quent generations. For homogeneity, we follow the nomenclature by 
Carretta et al. but we always consider all Na-rich stars as a single class 
(I+E), for simplicity. 
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NGC 


Alt. Name 


n„ 
arcmin 


foul 

arcmin 


rtn/r c 


rtn/n 


r BM /n, 


7* outl? t 


E(B-V) 


[Fe/H] 


2419 




0.0 


7.0 


0.0 


0.0 


7.3 


1.0 


0.11 


-2.12 


5024 


M53 


1.0 


7.0 


2.8 


0.9 


6.3 


0.3 


0.02 


-1.99 


5272 


M3 


0.0 


10.0 


0.0 


0.0 


8.9 


0.3 


0.01 


-1.57 


5466 




0.0 


10.0 


0.0 


0.0 


4.4 


0.3 


0.00 


-2.22 


5904 


M5 


0.0 


10.0 


0.0 


0.0 


4.7 


0.3 


0.03 


-1.34 


6205 


M 13 


1.0 


10.0 


1.3 


0.7 


6.7 


0.4 


0.02 


-1.54 


6341 


M92 


0.0 


10.0 


0.0 


0.0 


9.2 


0.7 


0.02 


-2.28 


7078 


M 15 


0.0 


10.0 


0.0 


0.0 


9.4 


0.5 


0.10 


-2.32 


7089 


M2 


1.0 


10.0 


2.9 


1.1 


10.7 


0.5 


0.06 


-1.62 



Notes. Global parameters a re from the most rece nt version of the lHarrisI dl996h on-line database (year 2003), e xcept for [Fe/H] f or both M 5 and 
M 15, which are taken from Carretta et al. (2009a), and structural parameters for NGC 2419 that are taken from Bellazzini ( 20071) . 



note that most of the RGB stars displaying the color segrega- 
tion as a function of Na abundance are fainter than the RGB 
bump (at g a 15.4, see Fig. 2, below), as expected if the ob- 
served chemical anomalies are not due to extra-mixing phenom- 
ena k nown to occur in evolutionary phases brigh ter than this fea- 
ture dGratton et al]l2000tlSmifh & Martellll2003l) . 

However, the clear u-g color segregation between Na-poor 
and Na-rich stars observed in M 5 is not seen in M 15 (lower 
panel of Fig. [TJ. A likely explanation of this different behavior 
calls into play the difference in overall metal content between 
the two clusters. The iron abundance in M 5 is ten times higher 
than in M 15 and the average abundance of the light elements 
should scale similarly. Hence, the same degree of N abundance 
with respect to iron (as expressed by [N/Fe]) corresponds to 
very different absolute abundances of N ([N/H]). This, in turn, 
should correspond to significant differences in the strength of 
the absorption features that are supposed to drive the spreads 
observed in broad -band near-ultraviolet photometry (see also 
Mart ell et al.ll200 8). In particular, the strengths of spectral lines 
of diatomic molecules (such as CN, NH) depend quadratically 
on the overall metallicity and it has been noted that CN bands 
become very weak in the spectra of GC giants for [Fe/H] 5= -1.8 
dSmifhll2002l) . In this context, we note that the clusters for which 
the color spread in the RGB has been detected using the broad U 
(or u, in the present case) filter and correlated with spectroscopic 
light-element abundances, have intermediate metallicities, i.e. 
[Fe/H]=-1.34 for M 5, [Fe/H]=-1.17 for M 4, [Fe/H]=-1.51 
for NGC 3201, rFe/H1=- 1.55 for NGC 6752 (metallicities 
from lCarrettaltal]|2009bl), and rFe/H l=-1.18 for NGC 1851, 
(metallicity from Carretta et al. 20 1 Pel) . 



On the other hand, we need also consider that the RGBs of 
M 5 and M 15 see m to display similar spreads in the (V,c v ) plane 
dYong et al.l2008l) . thus not supporting the above hypothesis. It is 
difficult to draw a firm conclusion at the present stage, in partic- 
ular if we consider that we still lack a detailed theoretical under- 
standing of the mechanism at the origin of the color spread. This 
kind of investigation is clearly beyond the scope of the present 
paper. For our present purposes, we take as the basis of the fur- 
ther analysis of the clusters in our sample the clear evidence that, 
at least in some cases, star-to-star differences in light-element 
abundances in the RGB of GCs can be discerned with the u-g 
color. 

In particular, the sought-after anomaly should be distinguish- 
able as a significant color spread in the RGB in the (g,u-g) CMD, 
not seen in the (g,g-r) CMD and unaccounted for by other ef- 
fects, such as an increase in the photometric error and/or an in- 



crease in the amplitude of the effect of differential extinction (see 
below). Independently of its true origin, we refer in the follow- 
ing section to the observational effect under consideration as a 
UV-spread (UV-s hereafter, for brevity) and we will look for its 
presence in the clusters of our sample, trying to establish, in the 
various cases, wether it is (at least partially) due to intrinsic star- 
to-star physical differences or it can be fully accounted for by 
trivial effects. 

3. Detecting anomalous u-g spreads 

In Fig. 2, we present the g,g-r CMDs, flanked by their g,u-g 
counterparts for all the clusters listed in Tab. [TJ In all cases, two 
horizontal segments enclose the magnitude range to which we 
limit our analysis of the color spread: we tried to select similar 
portions of the RGB in all clusters, possibly spanning a wide re- 
gion below the horizontal branch. Given the different distances 
and reddenings of the various clusters, this was not always pos- 
sible (in particular, for NGC 2419 which is far more distant than 
the other GCs in our sample). This gives rise to differences in the 
sensitivity of the adopted method as (a) the color spread corre- 
sponding to a given abundance spread is always observed to de- 
crease with increasing lu minosity along the RGB, virtually dis- 
appearing at the R GB tip dYong et alj|2008h iMarino et alJ l2008: 
iMilone et al.ll2010l) . and (b) the same portion of the RGB occurs 
at different apparent magnitudes in different clusters, which is 
indicative of different photometric accuracy. 

The line located approximately at the red edge of the RGB, 
within the two horizontal segments, is a ridge line following the 
curvature of the observed RGB, and is taken as a reference to 
compute the color spreads, that is the difference between the 
color of a given star and the color of the ridge line at the same 
magnitude, A fo /, where col = g - r or col - u- g. Both A g _ r and 
A„_ ? are computed only for stars selected as candidate RGB in 
the g,g - r CMDs (heavier dots in Fig. 2). To limit the effects 
of possible spurious u-g outliers, we considered only stars for 
which -0.2 < A„_g < 0.05; these limits were found to be appro- 
priate for all clusters except NGC 2419 for which we adopted 
-0.8 < A u - g < 0.05. 

There are several cases in Fig. 2 in which a conspicuous 
broadening of the RGB in u-g with respect to g-r is apparent. 
However, as anticipated above, a few factors unrelated to phys- 
ical differences among cluster stars may also (in principle) pro- 
duce this effect. These factors and the methods we adopted to 
keep their effect under control can be summarized as follows: 

1. Field contamination. The degree of contamination by 
Galactic fore/background stars is very modest, because 
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Fig. 2. g-r,g and g,u-g CMDs for the GCs in our sample. In some cases, the RGB is truncated due to SDSS saturation. The red 
arrows in the upper right corner of each panel are the reddening vectors whose amplitudes correspond to the average E(B-V) values 
reported in Table Q] Horizontal black arrows mark the position of the RGB bump. Stars selected as candidate RGB on the (g,g-r) 
CMDs are plotted as heavy dark dots, the remaining ones as lighter red dots. The curves approximately tracing the red edges of the 
RGBs are used as references to compute the color spread distributions shown in Fig. [3j the the two horizontal segments display the 
portion of the RGB that is used to compute those distributions. In most cases, note that the asymptotic giant branch starts above the 
bright end of the adopted selection box (NGC 2419 is an obvious exception). 



of the combination between the relatively high (absolute) 
Galactic latitude of the considered clusters (all have \b\ > 25° 
and four \b\ > 70°) and the relatively small area of the consid- 
ered annular fields. We used the Galactic model TRILEGAL 
dGirardi et alJl2005h to obtain a conservative estimate of the 
degree of contamination affecting the samples of candidate 
RGB stars considered in the present analysis (see Fig. 2). We 
found that the fraction of Galactic field stars in our samples 



is lower than 5% for seven of the nine clusters, reaching 8% 
for M 15 and 10% for M 92. Moreover, the effect of even 
this small degree of contamination should be minimized by 
the way in which we selected candidate RGB stars. In the 
following, we analyze the color spreads of candidate RGB 
stars selected in the g,g-r CMD as the most tightly clustered 
along the narrow RGB sequence in this plane (heavy dark 
dots in Fig. 2). This means that we compare the g-r and u- 
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g color spreads produced by exactly the same stars. Hence, 
any field star artificially broadening the distribution of A co / 
in u-g should also have a similar effect on g-r. 

2. Blendings and artifacts. Any photometry of crowded fields 
(especially ground-based seeing-limited ones) is expected to 
include some fraction of sources that are either blends of two 
(or more) fainter stars, or non-stellar sources, such as unre- 
solved distant galaxies or flukes in the halo of bright stars, 
both misclassified as stars. While it is likely that our samples 
of RGB candidates include some of these sources, their effect 
is expected to broaden the distribution of A eo ; by a similar 
amount in both of the considered colors, given the adopted 
selection, as for the case of contamination described above. 

3. Differential reddening. In principle, if there is a star-to-star 
difference in the degree of stellar extinction (due to spatial a 
variation in the reddening over the considered field of view) 
this will produce a larger color spread in u-g than in g-r, 
mimicking the effect we are looking for. However, the ratio 
of the expected spreads is quite small ^("j^ = 1.29, mean- 
ing that a difference of 0.02 mag in E(B-V) would corre- 
spond to 0.019 in E(g-r) and 0.024 in E(u-g), i.e. a mere 
difference of 5 millimag between g-r and u-g. The clusters 
considered here have, in general, low average reddening val- 
ues (four having E(B - V) < 0.02, seven E(B - V) < 0.06, 
and all nine E(B - V) < 0.11) and, as far as we know, no 
indication of differential reddening has ever been reported 
in the literature. In the cases for which we report the detec- 
tion of significant u-g spread, the full width at half maxi- 
mum (FWHM) of the distribution spread in u-g is 2 times 
larger than in g-r, more than the factor of 1 .29 that can be at- 
tributable to differential reddening alone. Finally, the case of 
M 5 presented here as well as the previous cases reported in 
the literature (see Sect. [TJ suggests that differential redden- 
ing does not play a major role in producing the UV-s, in the 
cases where this has been revealed up to now. We conclude 
that the effects of differential reddening for the considered 
sample, should be negligible. 

4. Photometric errors. Since all the SDSS observations are per- 
formed at fixed exposure time and RGB stars emit much less 
light in the near-UV than in the visible ranged, any given 
star in our sample has larger photometric errors in u than in 
either g or r: as a consequence, the color spread due to pho- 
tometric errors is larger in u-g than in g-r. We see below that 
this is the most serious problem in the present analysis, as 
there are cases in which a large observed UV-s can be fully 
accounted for by (relatively) large photometric errors that 
hide any (possible) underlying signal associated to real dif- 
ferences among stars. To take this effect into due account, we 
divided the color spread of each star by the associated pho- 
tometric uncertainty in each color A ml = ^ (normalized 
color spread), which should be (approximately) expressed 
in units of standard deviations. Comparing A , distributions 
in g-r and u-g, we can check whether there is any significant 
UV-s in addition to that due to photometric errors; if the lat- 



6 There are other factors concurring to lower the signal-to-noise ratio 
of u observations with respect to g or r ones. For example: (1) CCDs are 
less sensitive in near-UV than in visible bands; (2) as in the SDSS the 
photometry of a given field is taken simultaneously in all the ugriz pass- 
bands and the seeing worsens at shorter wavelengths, u images hence 
experiencing the worst seeing; (3) at fixed atmospheric transparency 
conditions and air mass, u light suffers from the highest amount of at- 
mospheric extinction. 



ter were the main contributor to the observed u-g spread, the 
two distributions should be indistinguishable. 

In Fig. [3] we plot the distribution of A co / and A co[ for all the 
considered clusters. Histograms in g-r are plotted as dotted lines, 
and those in u-g as continuous lines. The A' col distributions have 
been shifted by (small) arbitrary amounts to ensure that their 
maxima approximately coincide with A' , = 0, to allow a more 
direct comparison between the distributions in the two color in- 
dicefl We shall see below that the significance of the detected 
differences is maintained independently of the adopted shifts. 

It may be useful to start the discussion of Fig.[3]from the pair 
of panels in the lower right corner (NGC 2419). As readily vis- 
ible in the CMD of Fig. 2, there is an impressive broadening of 
the RGB color spread passing from g-r to u-g in this case. The 
broadening in u-g color is so large that we had to adopt a differ- 
ent horizontal scale in this panel to accommodate the bulk of the 
A„_ g distribution. However, the distributions of A co/ are indistin- 
guishable, with essentially all the stars lying within +3cr of the 
mean, as expected for (approximately) normal errorfl Hence, 
all the UV-s observed in NGC 2419 can be accounted for by 
the effect of photometric errors. We note that this does not mean 
that there is no underlying spread in light-element abundance, 
but that the photometric accuracy is not sufficient to reveal it, if 
it exists. The same is true for NGC 5466, while the case of M 15 
is discussed in more detail below. 

On the other hand, in all the cases considered in the left raster 
of panels pairs of Fig.[3](M 5, M 53, M 13, M 3, and M 2), plus 
M 92 in the right raster of panels, the differences in the A co ; 
distributions correspond to significant differences in the distri- 
butions of A ml , implying that the detected UV-s cannot be en- 
tirely accounted for by observational effects, thus requiring the 
presence of some physical star-to-star difference. In particular, 
the A„_ g distributions display extended tails toward the blue that 
are completely lacking in the remarkably symmetric g-r distribu- 
tions. In the case of M 5, we know from Sect.[2]that the observed 
u-g spread correlates with Na abundanc es (and it is likely to be 
caused by differences in N abundances. lYong et al.| [2008). It is 
interesting to note that the A co i distribut ion of this cluster appears 
to be bimodal (as in the case of M 4 Marin o et al . 2008), even if 
it is not possible to firmly establish the statistical significance of 
this feature. Hints of multi-modality are indeed visible in all the 
A co i (and A cg[ ) distributions of these six clusters. 

A Kolmogorov-Smirnov (KS) test would be the most 
straightforward non-parametric way to quantitatively establish 
the statistical significance of the detected differences between 
the u-g and g-r A col distributions. However, it is well known that 
this test is very sensitive to offsets between distributions and the 
shifts adopted in Fig. [3] are quite arbitrary. In principle, a given 
choice of the shift between the two distribution may add spuri- 
ous significance to their mutual difference measured by the KS 
test. To circumvent this problem, we proceeded as follows: (a) 
keeping the g-r distributions shown in Fig.[3]fixed, we moved the 
u-g distribution from -5 to +5 A col units in steps of 0. 1 ; (b) at 
each step, we performed the KS test computing the probability 
that the two samples are drawn from the same parent population, 
P%T', and, (c) we adopted the value of the shift that maximize 



7 In particular, the g-r distribution has been shifted to ensure that its 
well defined maximum occurs at A col = 0, while for the broader u- 
g distributions, with less clear peaks, we searched for a compromise 
between placing the peak at A col = and matching the right (red) edge 
of the two distributions, to make the comparison easier. 

8 we note that this is true for all the g-r distributions of A col . 
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Fig. 3. Distributions of u - g (solid histograms) and g-r (dashed histograms) color spreads (with respect to the RGB fiducials of 
Fig. 2) for the clusters in our sample. Panels marked with (a) show the distributions of absolute color spreads, while those marked 
with (b) show the normalized color spreads, in units of <x. 



Table 2. Dimensions of the samples and results of KS tests. 



NGC 


W stars 


Shift 


ppHot 
KS 


NuV-blue 


Nuv-red 


prad 
KS 


M5 


652 


0.2 


1.1 x lO" 11 


199 


453 


2.0 x 10-" 


M2 


492 


-0.3 


6.0 x 10" 10 


132 


360 


< 1.0 x 10-" 


M 3 


501 


1.5 


3.9 x 10-'° 


206 


295 


3.8 x lO" 10 


M 13 


596 


0.4 


5.7 x 10~ 8 


185 


411 


1.3 x 10~ 4 


M92 


442 


0.4 


3.8 x 10- 7 


139 


303 


< 1.0 x 10-" 


M53 


394 


0.4 


2.9 x 10~ 7 


75 


319 


5.7 x 10-" 


M 15 


371 





8.9 x 10~ 4 


110 


261 


7.8 x 10-" 


NGC 5466 


172 


0.2 


0.60 


37 


135 


0.15 


NGC 2419 


693 


-0.2 


0.79 


94 


599 


0.24 



Notes. N smrs : total number of candidate RGB stars selected for the analysis of color spread. "Shift" is the differential shift in A co( that is found to 



maximize P 



,phot 



Nrj 



and No 



are the number of stars having A' , < 2.0 and A' , > 2.0, respectively. P^T and P'£t are defined in the text. 



Clusters are listed in order of increasing P 



,phol 



Pg- , i.e. the shift that minimize the significance of the differ- 
ence between the two distributions. In this way, we are guaran- 
teed that the residual differences considered by the KS test are 
genuine differences in the shape of the distributions, not due to 
unphysical shifts. The adopted differential shifts and the corre- 
sponding values of P^' are reported in Tab. [2] The significance 
of the difference in the u-g and g-r distributions of A c ; is very 
high for all the clusters whose distributions are plotted in the 
left hand of Fig. [3]plus M 92 (P P K h ° r < 1(T 4 ). In particular, we 
note that the difference is also significant in the case of M 15: 
the probability that the u-g and g-r distributions of A' col mea- 
sured in this cluster are drawn from the same parent population is 
just 0.08%. A careful inspection of the A m[ distributions for this 
cluster in Fig. [3]reveals that while the wings virtually coincide, 
the core of the g-r distribution is far more peaked than its u-g 
counterpart, which exhibits hints of bi-modality (compare with 
the cases of NGC 5466 and NGC 2419, where the distributions 
nearly coincide in both the wings and the core). Hence, while 



the amplitude of the effect is insufficient to provide an obvious 
color segregation in the CMD of this cluster (see Fig.[TJ, the un- 
derlying signal is there and can be revealed once the much larger 
photometric sample is considered and the observational effects 
are properly taken into account, in agreement with the results of 
lYong et al] (120081) . This conclusion is strongly supported by the 
behavior of the radial distributions of RGB stars as a function 
of their color spread in this cluster, as discussed in the following 
section. 



We note that the four clusters displaying the most obvious 
and significant difference in their A' u _ g and A' g _ r distributions 
have metallicity around [Fe/H]= -1.5, while M 53, M 92, and 
M 15, which have larger P P KS ', are significantly more metal-poor 
([Fe/H]^ -2.0). This is likely due to the weakening of the intrin- 
sic UV-s effect at low metallicity, discussed in Sect. 12.11 
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3.1. The radial distribution of UV-red and UV-blue stars 

As anticipated in Sect. [T| it would be very interesting to check 
whether RGB stars with different u-g colors (at the same magni- 
tude) may have different radial distributions. According to some 
models of the early enrichment of GCs, at the end of the enrich- 
ment phase (lasting a few hundred Myr), stars from the first gen- 
eration (P) should be significantly less concentrated toward the 
center of the cluster than stars of later generatio n(s) (I+E), born 
from material polluted by the ashes of P stars dD'Ercole et al.l 
2008). While it may be expected that this difference should have 
been la rgely erased long ago by the dynamical e volution of the 
cluster dD'Ercole et al.l2008tlbecressin et al.l2008l) . the (sparse) 
available observational evidence about radial distributions of the 
two populations is in qualitative agree ment with the prediction of 
the above quoted model. For example, Kravtsov et al.l (I2.Q 1 Obllah 
found that the stars on the red side of of the broad U-I distribu- 
tion of the RGB of NGC 1261 and NGC 3201 are significantly 
more centrally conc entrated than those on the blue side. In the 
case of NGC 3201. ICarretta et al.l (l2010al) used detailed Na, O 
abundances from high resolution spectroscopy to show that Na- 
rich RGB stars are redder (in U-I) and more centrally concen- 
trat ed than N a-poor RGB stars. 

ICarretta et al.1 (|2()09b) merged the spectroscopic sample of 
RGB stars in 19 GCs, normalizing the distance of each star from 
the center of its cluster by the cluster r/,, and compared the ra- 
dial distributions of P, I, and E stars. They found that I stars are 
significantly more centrally concentrated that P stars, while E 
stars appear slightly l ess concentrated than P stars. However, as 
discussed in detail by ICarretta et al.l d2009b), this result may be 
affected by serious selection biases, inherent to the process of 
efficient fiber allocation in multi-object spectroscopy of Galactic 
GC stars. 

A difference in the radial distribution of the two sub giant 
branches of NGC 1851 ( whose stars are pr esumed to differ in 
terms of CNO abun dance; ICassisi etal . 2008) was detected by 
Zoccali et al. (2009), but the re sult was not fully con firmed by 
Milone et all d2009l) . However, ICarretta etal] (1201 Pel) find that 
the radial distribution of the RGB stars they analyzed depends on 
their iron abundance, more metal-poor stars being more centrally 
concentrated than their more metal-rich counterparts. 

Different populations with different r adial distributions are 
known to be present in oj Centauri (iPancino et al ] 120031: 
ISollima et ail 120071: IVillanova etal] 120071) . but this system is 
more complex and may have a dif ferent origin with from clas- 
sical GCs (see ICarretta et al.|l20 1 OcL and references therein). 

To follow-up this line of investigation with our sample, we 
divided the selected RGB candidates in each of the considered 
clusters into two sub-samples, according to the value of their 
normalized u-g spread A ,. For brevity, we dub UV-blue the 
stars having A , < -2.0 and UV-red those having A , > -2.0, 
in the scale of Fig. [3] (i.e. using the shifts adopted there). It is 
important to recall that any observational effect potentially af- 
fecting the radial distributions (notably, the radial variation in 
the degree of completeness, due to the increase in crowding to- 
ward the center) must affect the two sub-samples exactly in the 
same way, as they have the same distribution of magnitudes. 

The comparison between the radial distributions of UV-blue 
and UV-red RGB stars for all the considered clusters is presented 
in Fig. |4j the probability that the two distributions are drawn 
from the same parent population according to a KS test (P T £f) is 
reported in the last column of Tab. [2] In all the clusters for which 
we detected a significant intrinsic u-g spread (including M 15), 
the UV-red population is obviously more centrally concentrated 




0246B 02468 02468 
R [arcmin] R [arcmin] R [arcmin] 



Fig. 4. Comparison between the radial distribution of UV-blue 
(continuous blue line) and UV-red (dashed red line) for all the 
considered clusters. The data within the radius marked by the 
dotted lines in the plots of the distributions for M 53, M 2, and 
M 13 were not included in the analysis (see Sects. l2land 12. U . 
See Table [2] for the significance of the detected differences, as 
measured by the KS test. 



than the UV-blue one, with P™^ always lower (and in most cases 
much lower) than 0.02%. We note that this result is very insen- 
sitive to the actual choice of the A col threshold; the difference 
remains highly significant for a large range of adopted thresh- 
olds (in particular, for the whole range -3.0 < A col < 0.0). 

This result clearly provides further support to the physical 
significance of the u-g spread: it is very hard to conceive how 
any spurious observational effect can be associated with such a 
strong difference in the radial distribution. Moreover, it suggests 
that the higher degree of central concentration of UV-red (pu- 
tative I+E) stars with respect to UV-blue (putative P) stars may 
be a general characteristic of all GCs where intrinsic UV-s can 
be detected, therefore any model intended to explain the origin 
of the spread in light-element abundances in GCs should also be 
able to reproduce this feature. 

Fig. |4]aims to demonstrate the high level of statistical sig- 
nificance of the detected differences between the radial distribu- 
tions of UV-red and UV-blue stars. To allow a more direct com- 
parison with the pred i ctions of chemo-hydro-dyn amical models 
dD'Ercole et aDl2008l: IPecressin et al. 2008. 12010I) . in Fig. [5] we 
present the radial profile of the ratio of the number of UV-red to 
UV-blue stars ( WLzsi \ where the radial coordinate is expressed 
in units of the clusters core radii (r c ), and the distances of 1, 
2, 3, 4, and 5 half-light radii (r/,) are also indicated (vertical 
dashed lines). This ratio should approximately scale as the ratio 
of second generation(s) (I+E) to first generation (P) stars, whose 
radial distribution is a typical outcome of the considered mod- 
els. However, we stress that while the overall shape of the ob- 
served profiles does not change much when different thresholds 
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that specific model, after 25 half-light relaxation times of evolu- 
tion, the profile of the second-to-first generation number ratio is 
nearly flat within r 0.5 r/, (a region always entirely enclosed 
in the innermost bin of our profiles, or not even included in the 
considered sample, in the cases of M 2 and M 53), then has de- 
clined by a factor of a few by 2 r/,, at the limit of that figure. 
This is in fair agreement with the observed profiles^! a broad 
agreement is al s o foun d with some of the models presented in 
Decressin et al. see the middle bottom and right bottom 

panels of their Fig. 1). A detailed comparison with models is 
far beyond the scope of the present paper. On the other hand, 
we feel that Fig. [5] provides a very useful set of observational 
constraints that must be reproduced by models of GC formation. 
Unfortunately, our data do not sufficiently probe the innermost 
regions (r < 0.5 r/,) of the considered clusters; and complemen- 
tary HST observations are probably needed to check this part of 
the profiles. Finally, a fully meaningful comparison with the pro- 
files of Fig. [5] would require dedicated models for each cluster, 
taking into account its specific structural and dynamical proper- 
ties, as well as accounting for its evolution for a Hubble time. 



Fig. 5. Ratio of the number of UV-red to UV-blue stars as a func- 
tion of distance from the cluster center, for the seven clusters in 
which we detected a significant intrinsic UV-s. The radial co- 
ordinate is expressed in units of cluster core radii, the dashed 
vertical lines marks the radial distanc es corresp onding to 1,2, 
3, 4, 5 half-light radii (r c and r/, from[Harris 1996). The ratio is 
computed in radial bins 2.5' wide in steps of 0.5'. Note that the 
actual value of the ratio depends on the adopted threshold be- 
tween UV-red and UV-blue stars: here we adopted A co/ = -2.0 
as in Tab.|2]and Fig. [4] 

between UV-blue and UV-red stars are adoptecQ, the true values 
of Nuv-red/Nuv-biue depends on the adopted threshold, hence 
they cannot be directly compared with second-to-first genera- 
tion ratios computed elsewhere. The only safe conclusion that 
we can draw here is that UV-red stars are more abundant than 
UV-blue stars (in the innermost ~ 3 - 4 r/,) for any threshold 
A , < -1.0, albeit with cluster-to-cluster differences of a factor 
of a few; this is in qua litat ive agreement w i th the predictions by 
iD'Ercole et ail d2008l) andlDecressin et ail d2010h . and with the 
observations bv lCarretta et alJ d2010dl) . 

The profiles shown in Fig. [5] cover the range given by 
0.5 - 1 Oi an d > 4 r/,. In nearly all cases, the profiles display 
a relatively steep decline from the innermost bin out to a radius 
of ~ 3 Oi where they flatten at a level Nw-redlNuv-blue ~ 1, re- 
maining approximately flat out to the last observed point. Hints 
of another increase at r/, <: 10 r/, are seen in the profiles of M 3 
and M 2 but their significance seems only marginal, if any. A no- 
table exception to this general trend is provided by M 13, whose 
profile is nearly flat out to ~ 3 Oi, then declines into another flat 
branch at r ~ 4 r/,. 

It is interesting to compare the o bserved profile s with the 
model predictions shown in Fig. 18 of lD'Ercole et al. |(|200i).In 

9 Over the range -3.0 < A col < -1.0. In particular, the slope of the 
inner rising branch of the profile changes with the adopted threshold 
(the relative height of the central peak grows by a factor of S 3-4 
changing the threshold from -1 to -3) but the radius where the profile 
flattens remains unchanged. In Fig. [5] we have adopted A ( , o/ = -2, as 
above. 



3.2. Comments on individual clusters 

In this section, we briefly report on previous results about abun- 
dance and/or UV color spreads in the clusters considered in the 
present paper. 

Five of our c lusters were also considered in the study by 
lYong et all j2008h . namely M 5, M 3, M 13, M 15, and M 92. 
These clusters display a significant spread at all evolutionary 
stage s in the V, c y C MP (from photom etry by iGrundahl et al.l 
1 19991) . According to Yong et al. (2008), the observed spreads 
are comparable to that seen in NGC 6752, suggesting that all 
these clusters exhibit a ^2.0 dex dispersion in [N/Fe], i.e., the 
val ue found in NGC 6752. 

IGrundahl et alJ (1 19981) provided evidence of stars from 
two different populations in the horizontal branch of M 13. In 
particular, they noted that at any given (b - y)o color, there is a 
large spread in the cq index (defined as Co = c\ - 0.2E(b - y)) 
that they interpreted as being due to star-to-star-variations in 
CNO abundances. 

Turning to spectroscopic analyses, Mar tell et al.l d2008l) 
showed that M 53 has a broad but not strongly bimodal distribu- 
tion of CN band strength, with CN and CH band strengths that 
are anticorrelated for unevolved stars. 

For the well-studied cluster M 5, ISmith & Norrisl d!983l) 
showed that the cyanogen distribution is strongly bimodal 
for a sample of 29 stars near the RGB tip. The clas sic Na/O 
anti-correl ation has been foun d by Sne den et al.l d 19921) and con- 
firmed by llvans et al.l (|2001). This was further correlated with 
the CN strength index, 5S(3839): stars with larger CN indices 
also have larger Na and Al abundances and lower O abundances 
than stars with lower CN indices. Finally, iBrilev etlril d 19921) 
and ICohen et alJ d2002l) found strong (and anti-correlated) 
variations in the abundances of C and N of stars down to the 
base of the RGB ofM5. 



10 However, it is expected that the inner flat region of the profile would 
progressively extend as dynamical evolution proceeds beyond the 25 f r /, 
time lapse considered in that simulation (E. Vesperini, private commu- 
nication). 
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An extended Na/O anti-correl ation was found als o in the 
extremely metal-poor cluster M 92 (ISneden et al.ll 1 99 lb . as well 
as stro ng (and anticorrela ted) variations in the abundances of C 
and N (ICarbon et alJl982h . 

Spectroscopic observations revealed star-to-star variations 
in the abundances of the CNO gro up elements among the M 3 
giants dPilachowski & Snedenl200lh . with bot h oxygen-rich and 
oxygen-poor stars coexisting in the cluster dKraftet alj|1992t 
Cohen et al .ll2005h . 

LeeTfcoOO) claimed the existence of a C-N anticorrelation 
among sta rs on the lower RGB of M 15, although no bimodality 
is found (ICohen et al.1 l2005h . iKavser et all (120081) showed a 
weak indication of a bimodality in the CH-CN plane (two 
clumps separated at CN — 0.6), although the observational 
errors are large compared to the separation of the two clumps. 
Pancino et al. (2010, submitted to A&A) detected a clear and 
bimodal CN and CH anticor r elatio n among the unevolved 
stars measured bv lKavser et al.1 J2008). Finally, Smith & Mateo 
( fl990T) found a CN-CH band strength anticorrelation and a 
possible bi-modality in their sample of red giants in the cluster 
M2. 

We note that for all the clusters in which we detected a signif- 
icant intrinsic u-g spread, there was previous evidence of inho- 
mogeneity in their stellar population reported in literature, sup- 
porting our findings. 

Previous detections of differences in radial distributions be- 
tween different populations in GCs were limited to the cases of 
co Cen, NGC 1851, NGC 1261, and NGC 3201, as discussed in 
Sect. 13. II The results presented here more than doubles the num- 
ber of GCs where these differences are detected at a very high 
degree of significance. 

4. Summary & conclusions 

We ha ve used publicly available u,g,r photometry from I An et al] 
(2008) to search for anomalous spread in near UV color (u- 
g) along the RGB of nine high-latitude, low-reddening, and 
well populated Galactic GCs. This anomalous spread (UV-s) 
was detected before in some clusters, using colors including 
other broad/intermediate band filters, such as Johnson's U or 
Stromgren u, and it was shown to be associated with the well 
known spread in the abundan ce of light elements (C, N, O, Na, 
etc.; see lCarrett a et al. 2009b, and references therein). The main 
results of our analysis can be summarized as follows: 

1 . We have introduced a method to remove the effect of pho- 
tometric errors from the observed color spread, normalizing 
the color residual about a fiducial RGB line of each consid- 
ered star by the associated photometric error (A' cgl ). By com- 
paring the distribution of A u _ g with that of A g _ r of the same 
sample of candidate RGB stars, we revealed anomalous u- 
g spreads in seven of the nine clusters. The observed effect 
has a very high statistical significance in all seven cases: ac- 
cording to Kolmogorv-Smirnov tests, the probability that the 
observed A' and A r distribution are drawn from the same 
parent population are always lower that 0.1%, but typically 
much lower than this. The lack of detection of any significant 
intrinsic u-g spread in the two remaining clusters (NGC 2419 
and NGC 5466) may be due to a real lack of chemical spread 
but it may also be associated with insufficient photometric 
accuracy and/or radial sampling. NGC 2419 is far more dis- 
tant than any other cluster in the sample and its RGB is ob- 



served at much fainter apparent magnitudes, implying much 
larger photometric errors in the range of interest; for the 
same reason, the available photometry for this cluster sam- 
ples mostly the portion of the RGB between the HB and the 
tip, where the amplitude of the UV-s effect is known to de- 
crease. The radial range explored here for NGC 5466 is the 
smallest of the whole sample, r < 4.1 r&. There are indica- 
tions that the UV-s effect may be weaker at very low metal- 
licity and this may also hamper our ability to detect it in clus- 
ters with [Fe/H]^ -1.8. 

2. In the case of M 5 ([Fe/H]=-1.34), we have demonstrated 
that the Na abundance correlates with the u-g color along 
the RGB, in the same way as in other clusters studied in 
the literature: Na-rich cluster stars have redder u-g colors 
than Na-poor cluster stars of the same magnitude, while 
they are indistinguishable in g-r. The same effect is not seen 
in M 15 ([Fe/H]=-2.32): we attribute this to the weaken- 
ing of the UV-s with decreasing metallicity, which is proba- 
bly associated with the extr eme weakening of CN lines for 
[Fe/H]< -1.8 dSmithll2002l) . as CN and NH features in the 
region 3000A< A < 4000A are thought to be at the origin 
of the effect. However, the case of M 5 illustrates that light- 
element abundance variations can be traced with the u-g 
color as done before with U - B, U - V, U — I, and c v in- 
dices (see, for examplelYong et alJ 2008: Ma rino et alJl2008l ; 
lHan et alJl2009l;ICarretta et alj|2009bll2010dl) . 

3. Dividing the RGB stars of each cluster into UV-blue and UV- 
red subsamples, according to their A u _ f , we found that UV- 
red stars are more centrally concentrated than UV-blue stars 
in all the seven clusters in which we detected a significant 
intrinsic UV-s. Kolmogorv-Smirnov tests have proven that 
the difference in the radial distributions of the two groups are 
highly significant in all cases: the probability that UV-blue 
and UV-red stars be drawn from the same radial distribution 
is always lower than 0.02%. 

4. The radial profiles of the ratio of UV-red to UV-blue stars 
typically show an approximately linear decline from the first 
sampled point (at si 0.5 - 1 r/,) out to 3 - 4 r/,, where they 
flatten and remain approximately flat out to the last sampled 
point (at si 4-8 r/,). This behavior is in qualitative agreement 
with the predictions of most recent models of formation and 
chemical evolution of globular clusters. It is interesting to 
note that in M 13 the profile is flat in the innermost ~ 3 r/,, 
then declines at larger radii. The shape of the observed pro- 
file provides a quantitative basis to test these theoretical sce- 
narios, once specific models reproducing the present-day sta- 
tus of the considered clusters become available. 



The results presented here clearly suggest that the difference 
in the radial distribution of first and second generation stars may 
be a general characteristic of globular clusters. Moreover, it has 
been demonstrated that near UV photometry can be a very effi- 
cient tool to trace light element spreads in very large samples of 
RGB stars in clusters, and it is especially well suited to study- 
ing the radial distribution of the various cluster populations. It 
must be considered that the very encouraging results presented 
in this paper have been obtained with a relatively small (2.5 m 
aperture) ground-base d telescope, und er non-ideal seeing condi- 
tions (< 1.6" FWHM. IXn et alj|2008l) . A survey performed with 
larger telescopes and/or under good seeing conditions (e.g., with 
service mode observations) and complemented with HST obser- 
vations of the cluster cores would provide a completely new in- 
sight into the problem, and would surely contribute to shedding 
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light on the mysterious earliest phases of the evolution of globu- 
lar clusters. 
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